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TITLE OF THE INVENTION 

Simulator for a chemical mechanical polishing 
BACKGROUND OF THE INVENTION 
Field of the Invention 

5 The present invention relates to a simulator, and more particularly to a 

simulator for a chemical mechanical polishing (CMP) process used in the manufacture of 
semiconductor devices. 
Description of the Background Art 

The miniaturization and high integration of large scale integrated circuits (LSIs) 
10 creates the tendency that circuit patterns formed on the LSIs have a minimum pattern 
dimension of 0. 1 \i m. A circuit pattern of an LSI can be formed in a manner that writes 
a design circuit on a transfer mask for implementing it on a semiconductor substrate by 
laser or electron beam, and then performs a batch optical transfer of the transfer mask 
pattern onto the semiconductor substrate by a projection transfer apparatus. 
15 The resolution R of the transfer apparatus is given by the following expression: 

R = kl ^/NA 

where kl represents a process constant, X represents a waveform, and NA 
represents a numerical aperture. 

The circuit pattern is formed by optical transfer method as described above, and 
20 a transfer in a defocus state produces a blurred image, resulting in poor image forming 
performance. Here, the extent of focus to which a predetermined image forming 
performance can be maintained is referred to as "depth of focus (DOF)" and given by the 
following expression: 

DOF = k2X/NA 2 
25 where k2 represents a process factor. 
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In the present condition that the fabrication dimension approaches 0.1 \x m, the 
depth of focus that can be ensured optical theoretically is only about 0.3 [x m. 

On the other hand, repetitive processes such as selective etching and film 
formation are executed on the semiconductor substrate, and irregularities (substrate 
5 irregularities) occur on the surface of the semiconductor substrate. 

The occurrence of substrate irregularities was not a serious problem in such 
semiconductor devices that integration degree is low and substrate irregularities is smaller 
than the depth of focus. However, as the fabrication dimension is smaller, the substrate 
irregularities is recently larger than the depth of focus, making it difficult to obtain a 
10 predetermined image forming performance. 

The substrate irregularities can be eliminated by for example the following 
methods: one in which some dummy patterns irrelevant to a real circuit pattern are 
properly disposed to increase the bulk of lower portions (i.e., dummy pattern method); 
and another in which a semiconductor substrate is planarized by polishing so as to cut the 
15 irregularities generated thereon by chemical mechanical polishing (CMP). 

A general description of the planarization technique by CMP process is 
contained in, for example, "ULSI Lithography Technical Innovation," pp71-86, issued 
November 10, 1994 by Science Forum Corp. 

With the miniaturization and high integration of LSIs as stated above, a CMP 
20 process becomes increasingly critical. To effectively execute the CMP process, there is 
need for simulation in consideration of various parameters. However, heretofore there is 
no effective simulator. 
SUMMARY OF THE INVENTION 

It is an object of the present invention to provide a simulator that can execute 
25 simulations in consideration of various parameters in a CMP process. 
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The present invention is intended for a simulator for a chemical mechanical 
polishing process for planarizing a semiconductor substrate. The simulator receives a 
pattern density data containing information about a pattern density per unit region of a 
fabrication pattern in a pattern forming process of a semiconductor device, and first and 
5 second measured data about height distributions of irregularities on the semiconductor 
substrate that are measured before and after a chemical mechanical polishing process 
executed with respect to the pattern forming process. The first measured data is 
compared with a first calculated data about a two-dimensional distribution of 
irregularities on the semiconductor substrate before the chemical mechanical polishing 

10 process, which is calculated from the pattern density data. A least squares analysis is 
performed to obtain a first correlation coefficient, and a parameter fitting is performed 
such that square of the first correlation coefficient becomes a maximum. Also, the 
second measured data is compared with a second calculation data about a 
two-dimensional distribution of irregularities on the semiconductor substrate after the 

15 chemical mechanical polishing process, which is calculated from the pattern density data. 
A least squares analysis is performed to obtain a second correlation coefficient, and a 
parameter fitting is performed such that square of the second correlation coefficient 
becomes a maximum. 

One parameter fitting is accomplished by comparing the first calculated data 

20 about the two-dimensional distribution of irregularities on the semiconductor substrate 
before chemical mechanical polishing process with the first measured data about the 
height distribution of irregularities on the semiconductor substrate before chemical 
mechanical polishing process. Another parameter fitting is accomplished by comparing 
the second calculated data about the two-dimensional distribution of irregularities on the 

25 semiconductor substrate after chemical mechanical polishing process with the second 
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measured data after chemical mechanical polishing process. Therefore, the adjustment 
of parameters before chemical mechanical polishing process can be separated clearly 
from the adjustment of parameters after chemical mechanical polishing process, thus 
leading to the simulator that can consider various parameters. This provides the 
5 advantage that if process conditions is changed or a new apparatus is added, adjustment 
may be accomplished merely by making a fine adjustment of parameters. 

These and other objects, features, aspects and advantages of the present 
invention will become more apparent from the following detailed description of the 
present invention when taken in conjunction with the accompanying drawings. 
10 BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a flowchart to explain a method for simulating a CMP process 
according to a first preferred embodiment of the present invention; 

Fig. 2 is a block diagram to explain the configuration of a simulator for a CMP 
process in the first preferred embodiment; 
15 Figs. 3A to 3C are conceptual diagrams to explain a mesh adjusting 

operation; 

Fig. 4 is a flowchart to explain a mesh adjusting operation; 
Figs. 5 and 6 are diagrams illustrating schematically the relationship between 
pattern density and height of irregularities of a fabrication object surface; 
20 Fig. 7 is a flowchart to explain a method for simulating a CMP process 

according to a second preferred embodiment of the present invention; 

Fig. 8 a block diagram to explain the configuration of a simulator for a CMP 
process in the second preferred embodiment; 

Fig. 9 is a flowchart to explain a method for simulating a CMP process 
25 according to a third preferred embodiment of the present invention; 
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Fig. 10 is a block diagram to explain the configuration of a simulator for a CMP 
process in the third preferred embodiment; 

Fig. 1 1 is a flowchart to explain the operation of taking a CMP image; 

Fig. 12 is a diagram illustrating schematically the state that a polishing pad is 
5 pressed against a fabrication object surface before CMP; 

Fig. 13 is a flowchart to explain a method for simulating a CMP process 
according to a fourth preferred embodiment of the present invention; and 

Fig. 14 is a block diagram to explain the configuration of a simulator for a CMP 
process in the fourth preferred embodiment. 
10 DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Technical Idea of the Present Invention 

When manufacturing large scale integrated circuits (LSIs), a plurality of LSIs 
called "sub-chips," including a TEG (test element group) and a process monitor, are also 
formed on a semiconductor substrate in addition to a target LSI called "main chip." The 
15 plurality of sub-chips are formed in a region other than a region for forming the main 
chip. 

Not only a circuit pattern of the main chip but also circuit patterns of the 
plurality of sub-chips are formed in a transfer mask used in the individual processes for 
forming the main chip and sub-chips on the semiconductor substrate. 
20 To form the transfer mask, it is possible to use in the state that design data of 

the main chip and sub-chips are already synthesized. However, it is preferably split and 
processed because the actual LSI design data is tremendous. Therefore, the design data 
are often split in units of sub-chips and the split data are used sequentially. 

Splitting design data in units of sub-chips facilitates the handling of tremendous 
25 design data. The present inventor has reached a technical idea that the handling of 
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design data can be further facilitated by converting the design data in units of sub-chips to 
pattern density data of small scale. 

Preferred embodiments of a simulation method and simulator for a CMP 
process based on the above technical idea will be described below. 
5 In the following description, pattern density data is to be defined as follows. 

A pattern graphic is determined by design data of a semiconductor device. 
The pattern density is the percentage of the area of graphic components, i.e., pattern 
components, contained in a unit region of the pattern graphic. For example, if pattern 
components occupy one-half of the unit region that is a rectangular region of 100 //m 
10 square, the pattern density is 50 %. 

As used herein, the term "unit region" denotes a rectangular region 
corresponding to a one-mesh region for obtaining a pattern density two-dimensional 
distribution data. By an AND logic with the one-mesh region, an AND operation with 
the mesh region is executed so that the area of a pattern component is calculated to obtain 
15 its pattern density. Thus, the pattern density obtained per one-mesh region is pattern 
density data. 

To the whole area of an individual sub-chip, such calculation for each of the 
processes is made to obtain a pattern density data related to a respective process of the 
individual sub-chip. 

20 The following first to fourth preferred embodiments premise that a simulation 

of a CMP process is performed with the use of pattern density data per process, as 
described above. 

The simulator according to the present invention can be implemented by a 
computer system, and its software can be created with the use of algorithm of the 
25 simulation method according to the present invention. 
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A. First Preferred Embodiment 

A description will be made of a simulation method and simulator for a CMP 
process according to a first preferred embodiment of the present invention. 
A-l. Simulation Method and Simulator for CMP process 
5 A method for simulating a CMP process will be described by referring to the 

configuration of a simulator for a CMP process 1 shown in Fig. 2, and by using the 
flowchart shown in Fig. 1 . 

Referring to Fig. 2, the simulator for a CMP process 1 receives a pattern density 
data Dl per process from a pattern density data storage device 10, and receives a 
10 measured data D2 about height distributions before and after CMP per process from a 
height distribution measuring device 20. 

The height distribution measuring device 20 can be implemented by using, for 
example, an auto focus function of an exposure device used in the manufacture of 
semiconductor devices. 

15 That is, the exposure device has the function of irradiating obliquely laser light 

to a semiconductor substrate and observing its reflected light to measure the height of the 
substrate. The use of this function enables to take a two-dimensional distribution of the 
height of a fabricated pattern that is formed on a semiconductor substrate. 

Note that examples of the height distribution measuring device 20 should not be 

20 limited to the auto focus function of an exposure device. For example, an atomic force 
microscopy (AFM) may be used. 

As previously described, in the manufacturing processes of semiconductor 
devices, the processing such as selective etching and film formation is performed 
repetitively on a semiconductor substrate, so that in every process, irregularities occurs on 

25 the surface of the semiconductor substrate. One technique of eliminating such substrate 
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irregularities is CMP process. Therefore, a CMP process is executed every time one 
process is performed. 

The height distribution of irregularities on the semiconductor device differs 
before and after performing a CMP process. The height distribution measuring device 
5 20A measures a height distribution of irregularities on the semiconductor device before 
performing a CMP process and that after performing a CMP process, and provides the 
measured data to the simulator 1. 

In the simulator 1, based on a coordinate data contained in a pattern density data 
Dl, a pattern density two-dimensional distribution calculating part 111 expands the 
10 pattern density data such that a mesh data is arrayed in two dimensions to obtain a 
two-dimensional image. This provides a pattern density two-dimensional distribution 
image DP per process (step SI). 

The height distribution measuring device 20 provides a measured data D2 about 
the height distribution of irregularities on a semiconductor substrate before and after 
15 performing a CMP process. The measured data D2 is given as a two-dimensional 
distribution image of the condition of the irregularities on the semiconductor substrate. 
The mesh of the pattern density two-dimensional distribution image DP is not always 
identical with that of the height distribution measured data D2. In order to adjust such 
that the distribution image DP and measured data D2 have the same mesh, a mesh 
20 adjusting part 1 12 makes a mesh adjustment (step S2). 

The operation of mesh adjustment performed in the mesh adjusting part 112 
will be described by using Figs. 3A, 3B, 3C, and Fig. 4. 

Figs. 3A to 3C are diagrams illustrating schematically the processing in the 
mesh adjusting part 112. 
25 A Fourier image F in a Fourier space shown in Fig. 3A is subjected to a reverse 
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Fourier transform to obtain a reverse Fourier image R in a real space shown in Fig. 3B. 
If a certain mesh is added into the Fourier space as indicated by the broken line in Fig. 3A 
and then a value of zero is added to the image data on the added mesh, as shown in Fig. 
3C, the resultant reverse Fourier image R has a higher density than the image of Fig. 3B. 
5 For example, if new meshes are added on a Fourier space in order to increase 

the number of meshes to 2 m X2 n times meshes in two dimensions, a reverse Fourier 
transform to a real space causes an up sampling to 2 m X 2 n times meshes in two 
dimensions. The increasing rate of the number of meshes can usually be set to any 
desired value. It is preferable to set to a rate, such as 2 m X2 n times, at which it is 

10 possible to use FFT. 

The mesh adjusting part 112 executes the operation of mesh adjustment by 
using such mesh interpolation. 

More specifically, as shown in Fig. 4, the mesh adjusting part 112 performs a 
Fourier transform of the pattern density two-dimensional distribution image DP provided 

15 from the pattern density two-dimensional distribution calculating part 111 and the height 
distribution measured data D2 provided from the height distribution measuring device 20 
(step S211). 

This provides their respective Fourier images in the Fourier space, as described 
with reference to Fig. 3A. At this time, the respective number of meshes and size of 
20 meshes are found and the Fourier image having a lesser number of meshes is adjusted to 
that having a greater number of meshes. 

In general, the measured data D2 has a lesser number of meshes. Therefore, a 
mesh interpolation is performed such that the number of meshes in the measured data D2 
is adjusted to that in the pattern density two-dimensional distribution image DP. 
25 That is, new meshes are added around the Fourier image in the measured data 
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D2 (step S212). 

The values on the new meshes are then set to "0" in step S213. 

Subsequently, the Fourier image is subjected to a reverse Fourier transform 
(step S214), to obtain a reverse Fourier image, and the measured data D2 after the mesh 
5 adjustment and pattern density two-dimensional distribution image DP are reconstructed 
(step S2 15). 

Through the foregoing adjustment operation, the mesh of the pattern density 
two-dimensional distribution image DP matches the mesh of the height distribution 
measured data D2, and these two data can be compared to each other. 

10 Returning to Figs. 1 and 2, when the mesh adjustment completes in step S2, if 

the measured data after mesh adjustment is measured data before CMP process (i.e., 
before polishing), it is provided as a measured data D21 to a height distribution 
calculating part 113, together with the pattern density two-dimensional distribution image 
DP. If the measured data after mesh adjustment is measured data after the CMP process 

15 (after polishing), it is provided as a measured data D22 to a Fourier calculating part 114, 
together with the pattern density two-dimensional distribution image DP. 

The height distribution calculating part 113 calculates a height distribution 
based on the pattern density two-dimensional distribution image DP, to obtain a height 
distribution data DPI about a fabrication object surface before CMP process (step S3). 

20 Here, a method for calculating a height distribution based on a pattern density 

two-dimensional distribution image DP will be schematically explained by using Figs. 5 
and 6. 

Fig. 5 shows a state that already fabricated circuit pattern PT and plurality of 
circuit patterns PT2 are arrayed on a semiconductor substrate SB. In Fig. 5, region Rl 
25 represents the region where the circuit pattern PI is disposed; region R2 represents the 
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region where the plurality of circuit patterns PT2 are disposed; and region R3 represents 
the region where no circuit patter is disposed. 

In region Rl, the circuit pattern PT1 is formed so as to cover the entire region 
and its pattern density is 1 00 %. In region R2, a 50 % of its entire region is covered with 
5 the circuit patterns PT2, and its pattern density is 50 %. The pattern density of region 
R3 is 0 %. 

In the manufacturing processes of LSIs, the process of forming an insulating 
film and metal film and the process of patterning these films are repetitively performed. 
Therefore, the insulating film or metal film (hereinafter referred to as a "laminated film") 

10 is to be formed on an already fabricated circuit pattern. When forming a laminated film, 
the material of the laminated film per unit area is supplied nearly uniformly throughout 
the entire semiconductor substrate surface. 

Fig. 6 shows a state that a laminated film SFM is formed on the semiconductor 
substrate SB, and there occurs a height distribution in the laminated film SFM due to a 

15 difference in pattern density between the circuit patterns. 

Let, dl is the formation thickness of the laminated film SFM, and d2 is the 
thickness of the circuit pattern PT1 or PT2. In region Rl having a pattern density of 
100 %, the total height HI of the laminated film SFM and circuit pattern PT1 is given by 
the following expression: HI = dl + d2 X 1 .0. In region R2 having a pattern density 

20 of 50 %, reflow and annealing processing fill grooves to facilitate planarization (the 
planarization effect during the film formation), so that the total height H2 of the laminated 
film SFM and circuit pattern PT2 is given by the following expression: H2 = dl + d2 
X 0.5. In region R3 having a pattern density of 0 %, the total height H3 of the 
laminated film SFM and circuit pattern PT2 is given by the following expression: H3 = 

25 dl + d2 X 0. In the above three expressions, the factors "1.0", "0.5", and "0" are 
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pattern densities. 

Since the formation thickness dl of the laminated film SFM is common to all 
the regions, it is relatively meaningless and can be eliminated. Accordingly, the height 
of each region can be given by the expression: d2 X (pattern density). 
5 The thickness d2 of the circuit pattern PT1 or PT2 is a process parameter that 

varies depending on the type of the pattern. 

Therefore, the height distribution calculating part 113 can obtain a height 
distribution of the fabrication object surface before CMP process in such a simple 
arithmetic of multiplying the pattern density two-dimensional distribution image DP by 
10 the thickness of a laminated film to be formed in the following next process. 

Returning to Figs. 1 and 2, first, the operation of parameter fitting before CMP 
process will be described. 

After taking the height distribution data DPI about the fabrication object 
surface before CMP process in step S3, the measured data before CMP process D21 and 
15 height distribution data DPI are provided to a correlation coefficient calculating part 1 1 8. 

In the correlation coefficient calculating part 1 18, the measured data D21 and 
height distribution data DPI are subjected to a least squares analysis to calculate a 
correlation coefficient (step S4). 

As used herein, the term "least squares analysis" refers to the following 
20 technique that two height distribution data are compared to each other and their similarity 
is analyzed by least squares method. 

Following is a brief description of the lease squares method. In data x and 
data y, their respective sample variances are expressed by the following equation (1) and 
equation (2), respectively. 
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S 2 = £ .-(I) 
x n 



S 2 .^-?)' -..(2) 



The sample covariance can be expressed by the following equation (3): 



Sxy = *(x-x)(y-y) ... (3) 
n 



In the above equations (1) to (3), x and y represent a mean value of data x 
10 and data y, respectively, and n represents the number of data. 

A correlation coefficient r given by the following equation (4) can be defined 
from the above-mentioned sample variances and sample covariance. 

r = — (4) 

SxSy 

15 The correlation coefficient calculated in step S4 has a value from -1 to 1. For 

a positive correlation, "1" represents a complete match. For a negative correlation, 
represents a complete match. Accordingly, it can be said that in the positive and 
negative correlations, a larger square value increases the degree of match between data x 
and data y, namely between the height distribution in the measured data X21 and the 

20 height distribution in the distribution data DPI. The maximum square of correlation 
coefficient indicates that the former height distribution is approximately the same as the 



14 

latter. 

Using the value of the square of correlation coefficient as an index, a parameter 
fitting part 1 19 performs a parameter fitting such that the square of correlation coefficient 
approaches the maximum (step S 1 0). 
5 Concretely, when forming a film on a semiconductor substrate (a fabrication 

object surface) having a certain pattern, the thickness of the film formed is a fitting 
parameter. There is of course other fitting parameter than the thickness of a laminated 
film. 

Here, a height distribution after forming the above film on the semiconductor 
10 substrate is measured by the height distribution measuring device 20, and the measured 
result is measured data D21. A calculated value based on a pattern density data when 
forming the above film is height distribution data DPI about the fabrication object surface 
before CMP process. 

Therefore, the formation thickness d2 of the laminated film SFM that was set in 
15 step S3 is changed so as to approach the measured data D21, i.e., increase the correlation 
coefficient. This is one of the parameter fittings before CMP process. 

Following is the operation of parameter fitting after CMP process. 
The pattern density two-dimensional distribution image DP that is provided to 
the Fourier calculating part 114 is then subjected to a Fourier transform, to obtain a 
20 two-dimensional Fourier image. This Fourier transform causes a projection from a real 
space to a frequency space, so that a two-dimensional image in the real space is 
transformed to a two-dimensional Fourier image represented by the magnitude of space 
frequencies (step S5). 

Here, a component having a small space frequency corresponds to a region in 
25 the real space where many isolated patterns are present, and a component having a large 
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space frequency corresponds to a region in the real space where many dense patterns are 
present. 

Subsequently, in a spatial filter part 115, the two-dimensional Fourier image is 
subjected to a spatial filter that permits only passage of components having a small space 
5 frequency. As the result, the components of small space frequency are selected and 
components of large space frequency are removed (step S6). The technique of spatial 
filter is well known. 

Here, the component having a small space frequency corresponds to a 

component that is a factor contributing to the phenomenon having a long correlation 
10 distance. The component having a long space frequency corresponds to a component 

that is a factor contributing to the phenomenon having a short correlation distance. 

Therefore, the spatial filter removes the components of large space frequency, 

leaving the components of small space frequency, i.e., the components that is the factor 

contributing to the phenomenon having a long correlation distance. 
15 When patterns of the same size exist in different densities in a CMP process, 

there occurs such a phenomenon that its polish velocity differs depending on the location. 

This phenomenon has an extremely long correlation distance, as long as 10 m to 100 

ju m, in some cases. 

Subsequently, a reverse Fourier calculating part 116 performs a reverse Fourier 
20 transform of the two-dimensional Fourier image holding only the components of small 
spatial frequency, thereby obtaining a reverse Fourier image, i.e., a pattern density 
two-dimensional distribution image DPX in the real space (step S7). 

This pattern density two-dimensional distribution image DPX indicates only the 
components that are the factor contributing to the phenomenon having a long correlation 
25 distance. This is a two-dimensional distribution image suitable for analyzing the 
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phenomenon having a long correlation distance. 

The above-mentioned two-dimensional distribution image DPX and measured 
data after CMP process D22 are then provided to the height distribution calculating part 
117. Based on the pattern density two-dimensional distribution image DPX, the height 
5 distribution calculating part 117 calculates a height distribution to obtain a height 
distribution data DP2 containing only factors that can cause the phenomenon having a 
long correlation distance (step S8). 

This height distribution calculation method need not be described herein 
because it is the same as the method for calculating the height distribution based on the 
10 pattern density two-dimensional distribution image DP, which is described previously 
with reference to Figs. 5 and 6. 

After the height distribution data DP2 is obtained in step S8, the data DP2 and 
measured data after CMP process D22 are provided to a correlation coefficient 
calculating part 118. 

15 The correlation coefficient calculating part 118 performs a least squares 

analysis of the measured data D22 and height distribution data DP2 to calculate a 
correlation coefficient (step S9). The operation of step S9 need not be described herein 
because it is the same as the analysis operation in step S4. 

Using the correlation coefficient obtained in step S9, as an index, the parameter 

20 fitting part 119 performs a parameter fitting such that the square of the correlation 
coefficient approaches the maximum (step S10). 

Here, assuming that a film is formed on a semiconductor substrate (a fabrication 
object surface) having a certain pattern, the height distribution measuring device 20 
measures a height distribution at the stage where this film is already polished by CMP. 

25 The measured result is measured data D22. On the other hand, a height distribution 
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obtained based on the two-dimensional distribution image indicating only the components 
that can cause the phenomenon having a long correlation distance is height distribution 
dataDP2. 

Therefore, the formation thickness d2 of the laminated film SFM that was set in 
5 step S8 is changed so as to approach the measured data D22, i.e., increase the correlation 
coefficient. This is one of the parameter fittings after CMP process. 

For example, the use of a two-dimensional distribution data after being 
subjected to a Fourier analysis enables to consider elastic deformation before a polishing 
is started with a polishing pad of CMP pressed against irregularities of a semiconductor 
10 substrate. This permits an analysis free from any influence on parameters such as 
polishing time and the number of revolutions of the polishing pad. 

The foregoing operations in steps SI to S10 are repeated with respect to 
measured data before and after CMP process in all the inputted manufacturing processes. 
A-2. Effects 

15 According to the method for simulating a CMP process and the simulator for a 

CMP process in the first preferred embodiment, before and after CMP process, a 
measured data and a simulation data are compared to obtain a correlation. Therefore, 
the adjustment of parameters before CMP process can be separated clearly from the 
adjustment of parameters after the CMP process, thus leading to a simulator that can 

20 consider various parameters. 

This provides the characteristic feature that if process conditions is changed or a 
new apparatus is added, adjustment may be accomplished merely by making a fine 
adjustment of parameters. 

It is also possible to know the degree to which the elastic deformation of a 

25 polishing pad affects the measured data D22, by monitoring the correlation between the 
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measured data after CMP process D22 and the height distribution data DP2 obtained 

based on the two-dimensional distribution image indicating only factors that can cause the 

phenomenon having a long correlation distance. 

B. Second Preferred Embodiment 
5 Following is a method for simulating a CMP process and a simulator for a CMP 

process according to a second preferred embodiment of the present invention. Similar 

reference numerals have been used in the following figures to denote similar parts that 

have the same configuration as in Figs. 1 and 2, and need not be described herein. 

B-l. Simulation Method and Simulator for CMP Process 
10 A method for simulating a CMP process will be described by referring to the 

configuration of a simulator for a CMP process 2 shown in Fig. 8, and by using the 

flowchart shown in Fig. 7. 

Referring to Fig. 8 5 the simulator 2 receives a pattern density data Dl per 

process from a pattern density data storage device 10, and receives a measured data D3 
15 about a height distribution of an under laid layer and height distributions before and after 

CMP per process from a height distribution measuring device 20. 

As used herein, the term "height distribution of the under laid layer" denotes the 

height distribution of the under laid layer of a fabrication pattern to be formed. 

Specifically, in the manufacturing processes of semiconductor devices, some processes 
20 such as selective etching and film formation are repeatedly performed with the use of 

about 30 types of masks. Accordingly, the semiconductor substrate surface before being 

subjected to a patterning can be assumed to be flat only at the time of its initial process. 

In the succeeding processes, some irregularities are already present on the semiconductor 

substrate surface before performing a patterning. Such irregularities on the 
25 semiconductor substrate is defined as a height distribution of an under laid layer. 
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In the simulator 2, based on a coordinate data contained in the pattern density 
data Dl, a pattern density two-dimensional distribution calculating part 211 expands the 
pattern density data such that a mesh data is arrayed in two dimensions to obtain a 
two-dimensional image. This provides a two-dimensional distribution image DP per 
5 process (step S21). 

In the measured data D3 about the height distribution of irregularities on the 
under laid layer and the height distributions of irregularities on the semiconductor device 
before and after performing a CMP process, which are provided from the height 
distribution measuring device 20, the condition of the irregularities on the semiconductor 

10 device is provided as a two-dimensional distribution image. However, the mesh of the 
pattern density two-dimensional distribution image DP is not always identical with that of 
the height distribution measured data D3. In order to adjust such that the 
two-dimensional distribution image DP and measured data D3 have the same mesh, a 
mesh adjusting part 212 makes a mesh adjustment (step S22). 

15 Through the foregoing adjustment operation, the mesh of the pattern density 

two-dimensional distribution image DP matches the mesh of the height distribution 
measured data D3, and these two data can be compared to each other. 

The mesh adjustment operation in the mesh adjusting part 212 is the same as 
that of the mesh adjusting part 112 shown in Fig. 2, and need not be described herein. 

20 At the completion of the mesh adjustment in step S22, the resulting measured 

data before it is subjected to a CMP process (i.e., before polishing) is provided as a 
measured data D31 to a height distribution calculating part 213, together with the pattern 
density two-dimensional distribution image DP. 

The measured data after CMP process (after polishing) is provided as a 

25 measured data D32 to a Fourier calculating part 215, together with the pattern density 
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two-dimensional distribution image DP. 

The height distribution calculating part 213 calculates a height distribution 
based on the pattern density two-dimensional distribution image DP, to obtain a height 
distribution data DPI about a fabrication object surface before CMP process (step S23). 
5 The method for calculating the height distribution based on the pattern density 
two-dimensional distribution image DP is described previously with reference to Figs. 5 
and 6, and need not be described herein. 

A measured data adding part 214 adds data of the height distribution of the 
under laid layer to the height distribution data DPI about the fabrication object surface 
10 before CMP process, to obtain a height distribution data DPI 1 (step S24). 

Following is the operation of parameter fitting before CMP process. 

After step S24 acquires the height distribution data DPI 1 about the fabrication 
object surface before CMP process that contains the height distribution data about the 
under laid layer, the measured data before CMP process D3 1 and height distribution data 
15 DPI 1 are provided to a correlation coefficient calculating part 220. 

The correlation coefficient calculating part 220 performs a least squares 
analysis of the measured data D21 and height distribution data DPI, to calculate a 
correlation coefficient (step S25). The analysis operation of step S25 is the same as that 
of step S4 shown in Fig. 1, and need not to be described herein. 
20 Here, assuming that a film is formed on a semiconductor substrate (a fabrication 

object surface) having a certain pattern, the thickness of the formed film is a fitting 
parameter. The height distribution after forming this film on the semiconductor 
substrate is measured by the height distribution measuring device 20, and the measured 
result is measured data D3 1 . The value calculated based on a pattern density data when 
25 forming this film is height distribution data DPI 1 . 
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Therefore, the formation thickness d2 of a laminated film SFM that was set in 
step S23 is changed so as to approach the measured data D31, i.e., increase the correlation 
coefficient. This is one of the parameter fittings before CMP process. 

Following is the operation of parameter fitting after CMP process. 
5 The pattern density two-dimensional distribution image DP provided to the 

Fourier calculating part 215 is subjected to a Fourier transform (step S26) and then 
subjected to a spatial filter in a spatial filter part 216 (step S27). The resulting image is 
then subjected to a reverse Fourier transform in a reverse Fourier calculating part 217, to 
obtain a reverse Fourier image, i.e., a pattern density two-dimensional distribution image 
10 DPX in the real space (step S28). 

The above-mentioned two-dimensional distribution image DPX and measured 
data after CMP process D32 are then provided to the height distribution calculating part 
218. Based on the pattern density two-dimensional distribution image DPX, the 
calculation part 218 calculates a height distribution to obtain a height distribution data 
15 DP2 containing only factors that can cause the phenomenon having a long correlation 
distance (step S29). 

The method for calculating the height distribution based on the pattern density 
two-dimensional distribution image DPX need not be described here because it is the 
same as the method for calculating the height distribution based on the pattern density 
20 two-dimensional distribution image DP, which is described previously with reference to 
Figs. 5 and 6. 

A measured data adding part 219 adds data of the height distribution of an 
under laid layer to the pattern density two-dimensional distribution image DPX, to obtain 
a height distribution data DP21 containing data of the height distribution of the under laid 
25 layer (step S30). 
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After step S30 acquires the height distribution data DP21 5 the data DP21 and 
measured data after CMP process D32 are provided to a correlation coefficient 
calculating part 220. 

The correlation coefficient calculating part 220 performs a least squares 
5 analysis of the measured data D32 and height distribution data DP21, to calculate a 
correlation coefficient (step S31). The analysis operation of step S31 is the same as that 
of step S25. 

Using the correlation coefficient obtained in step S3 1, as an index, a parameter 
fitting part 221 performs a parameter fitting such that the square of the correlation 

10 coefficient approaches the maximum (step S32). 

Here, assuming that a film is formed on a semiconductor substrate (a fabrication 
object surface) having a certain pattern, the height distribution measuring device 20 
measures a height distribution at the stage where this film is already polished by CMP. 
The measured result is measured data D32. On the other hand, a height distribution 

15 obtained based on the two-dimensional distribution image indicating only factors that can 
cause the phenomenon having a long correlation distance is height distribution data DP21. 

Therefore, the formation thickness d2 of a laminated film SFM that was set in 
step S29 is changed so as to approach the measured data D32, i.e., increase the correlation 
coefficient. This is one of the parameter fittings after CMP process. 

20 The foregoing operations in steps S21 to S32 are repeated with respect to 

measured data before and after CMP process in all the inputted manufacturing processes. 
B-2. Effects 

According to the method for simulating a CMP process and the simulator for a 
CMP process in the second preferred embodiment, before and after CMP process, a 
25 measured data and a simulation data are compared to obtain a correlation. Accordingly, 
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the adjustment of parameters before CMP process can be separated clearly from the 
adjustment of parameters after CMP process, thus leading to the simulator that can 
consider various parameters. 

This provides the characteristic feature that if process conditions is changed or a 
5 new apparatus is added, adjustment may be accomplished merely by making a fine 
adjustment of parameters. 

In addition, with the configuration that the measured data of the height 
distribution of the under laid layer in the individual process is added to the pattern density 
data, it is possible to consider the influence of the previous process and therefore permit a 
10 simulation suitable for manufacturing semiconductor devices having a laminated 
structure. 

C. Third Preferred Embodiment 

Following is a method for simulating a CMP process and a simulator for a CMP 
process according to a third preferred embodiment of the present invention. Similar 
15 reference numerals have been used in the following figures to denote similar parts that 
have the same configuration as in Figs. 1 and 2, and need not be described herein. 
C-l. Simulation Method and Simulator for CMP Process 

A method for simulating a CMP process will be described by referring to the 
configuration of a simulator for a CMP process 3 shown in Fig. 10, and by using the 
20 flowchart shown in Fig. 9. 

Referring to Fig. 10, the simulator 3 receives a pattern density data Dl per 
process from a pattern density data storage device 10, and receives a measured data D2 
about height distributions before and after CMP per process from a height distribution 
measuring device 20. 

25 In the simulator 3, based on a coordinate data contained in the pattern density 
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data Dl, a pattern density two-dimensional distribution calculating part 311 expands the 
pattern density data such that a mesh data is arrayed in two dimensions to obtain a 
two-dimensional image. This provides a two-dimensional distribution image DP per 
process (step S41). 

5 In the measured data D2 about the height distributions of irregularities on an 

under laid layer before and after performing a CMP process, which are provided from the 
height distribution measuring device 20, the condition of the irregularities on the 
semiconductor device is provided as a two-dimensional distribution image. However, 
the mesh of the pattern density two-dimensional distribution image DP is not always 

10 identical with that of the height distribution measured data D2. In order to adjust such 
that the two-dimensional distribution image DP and measured data D2 have the same 
mesh, a mesh adjusting part 312 makes a mesh adjustment (step S42). 

The mesh adjustment operation in the mesh adjusting part 312 is the same as 
that of the mesh adjusting part 1 12 shown in Fig. 2, and need not be described herein. 

15 Through the foregoing mesh adjustment operation, the mesh of the pattern 

density two-dimensional distribution image DP matches the mesh of the height 
distribution measured data D2, and these two data can be compared to each other. 

At the completion of the mesh adjustment in step S42, the' resulting measured 
data D21 before it is subjected to a CMP process (before polishing) and measured data 

20 after CMP process (after polishing) D22 are provided to a height distribution calculating 
part 313, together with the pattern density two-dimensional distribution image DP. 

The height distribution calculating part 313 calculates a height distribution 
based on the pattern density two-dimensional distribution image DP, to obtain a height 
distribution data DPI about a fabrication object surface before CMP process (step S43). 

25 The method for calculating the height distribution based on the pattern density 
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two-dimensional distribution image DP is described previously with reference to Figs. 5 
and 6, and need not be described herein. 

The height distribution data DPI about the fabrication object surface before 
CMP process and measured data before CMP process D21 are provided to a correlation 
5 coefficient calculating part 315, and these data and the measured data after CMP process 
D22 are provided to a CMP image calculating part 314. 

Following is the operation of parameter fitting before CMP process. 

The correlation coefficient calculating part 315 performs a least squares 
analysis of the measured data D21 and height distribution data DPI, to calculate a 
10 correlation coefficient (step S44). The analysis operation of step S44 is the same as that 
of step S4 shown in Fig. 1, and need not to be described herein. 

Subsequently in the parameter fitting part 316, the formation thickness d2 of a 
laminated film SFM that was set in step S43 is changed so as to approach the measured 
data D22, i.e., increase the correlation coefficient (step S47). This is one of the 
15 parameter fittings after CMP process. 

Following is the operation of parameter fitting after CMP process. 

The CMP image calculating part 314 calculates a two-dimensional distribution 
data after polishing that is obtained by CMP process, namely a CMP image, from the 
height distribution data DPI about the fabrication object surface before CMP process, as 
20 well as mechanical parameters such as Young's modulus and elastic coefficients of a 
polishing pad used in the CMP process (step S45). 

The operation of taking a CMP image will be described with reference to Figs. 
11 and 12. 

Referring to Fig. 11, first, in step S451 the CMP image calculating part 314 
25 calculates the shape of a polishing pad PAD when it is pressed against a fabrication object 
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surface, based on the height distribution data DPI about the fabrication object surface 
before CMP process. 

Fig. 12 illustrates schematically the state that the polishing pad PAD is pressed 
against the fabrication object surface before CMP process. Similar reference numerals 
5 have been used in Fig. 12 to denote similar parts that have the same configuration as in 
Fig. 6, and need not be described herein. 

Referring to Fig. 12, when the polishing pad PAD is pressed against the 
fabrication object surface, a characteristic phenomenon occurs at the boundary portions 
between regions Rl and R2 and that between regions R2 and R3. Specifically, the 
10 polishing pad PAD is brought into contact with the laminated film SFM at locations 
indicated by character "A" in these boundary portions, so that these locations are 
subjected to large stress and the laminated film SFM is well polished. On the other hand, 
at locations indicated by character "C", the polishing pad PAD is away from the 
laminated film SFM and these locations are subjected to less stress, resulting in poor 
15 polishing to the laminated film SFM. 

The height distribution of the laminated film SFM before polishing, i.e., the 
shape of the laminated film SFM, is given by the product of a pattern density 
two-dimensional distribution and thickness d2 of a circuit pattern PT1 or PT2. The 
two-dimensional distribution image of the polishing pad PAD is obtainable by the product 
20 of a reverse Fourier image and thickness d2. 

Returning to Fig. 11, after calculating the two-dimensional distribution image 
of the irregularities of the polishing pad (i.e., the pad shape), a two-dimensional 
distribution of stress exerted on the polishing pad is calculated (step S452). 

The stress exerted on the polishing pad PAD will be described with reference to 

25 Fig. 12. 
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Referring to Fig. 12, when the polishing pad PAD is pressed against the 
fabrication object surface, the polishing pad PAD is brought into contact with the 
laminated film SFM at locations indicated by character "A" in these boundary portions, 
whereas at locations indicated by character "C" are subjected to less stress, resulting in 
5 poor polishing to the laminated film SFM. 

Specifically, since large stress is being exerted on the regions, such as regions 
Rl and R2, where the change amount of the polishing pad PAD is large, these regions are 
polished promptly. In contrast, in the region free from any distortion such as region R3, 
and in the case that distortion occurs in the opposite directions in the region between 
10 projected patterns, the stress is zero, making it difficult to polish these regions. Even 
projected patterns of the same size will be subjected to different stresses and different 
polish velocities, depending on whether there is any other projected pattern that supports 
adjacent to these projected patterns. 

In step S452, such stress that varies depending on the location of the polishing 
15 pad PAD is calculated to obtain a two-dimensional distribution image of the stress. 

The two-dimensional distribution image of the stress exerted on the polishing 
pad PAD can be found from a difference between the shape of the laminated film SFM 
and the shape of the polishing pad PAD (i.e., the two-dimensional image of the 
irregularities). 

20 That is, the two-dimensional distribution image of the stress exerted on the 

polishing pad PAD is obtainable by multiplying a numeric value, which is obtained by 
subtracting the numeric data of the shape of the polishing pad PAD from the numeric data 
of the shape of the laminated film SFM, by Young's modulus (elastic coefficients). 

In step S453, a two-dimensional data DP3 about the irregularities on the 

25 fabrication object surface after polishing is calculated based on the two-dimensional 
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distribution image of the stress on the polishing pad PAD. 

The height of the fabrication object surface after polishing can be obtained as 
follows. First, a polishing amount (A) of a calculation object area is obtained by 
multiplying the following items: (i) a polishing rate (A/sec) that is determined by the 
5 material of the fabrication object surface, the material of the polishing pad and the 
number of revolutions of the polishing pad, etc.; (ii) a stress value (pascal) exerted on the 
polishing pad in the calculation object area; and (iii) a polishing time (sec). Next, the 
obtained polishing amount is subtracted from the height of the calculation object area of 
the fabrication object surface before polishing. 
10 Since in step S453 the area subjected to the highest stress, namely, the 

uppermost projecting area on the fabrication object surface, is large in polishing rate, it is 
possible to express the situation that the substrate is being planarized. This provides a 
two-dimensional distribution data DP3 of the irregularities on the fabrication object 
surface after polishing. 

15 Returning to Figs. 9 and 10, the two-dimensional distribution data DP3 and 

measured data after CMP process D22 are provided to the correlation coefficient 

calculating part 315. 

The correlation coefficient calculating part 315 performs a least squares 

analysis of the measured data D22 and two-dimensional data DP3, to calculate a 
20 correlation coefficient (step S46). The analysis operation of step S46 is the same as that 

of step S4 described with reference to Fig. 1, and need not be described herein. 

Using the correlation coefficient obtained in step S46, as an index, the 

parameter fitting part 316 performs a parameter fitting such that the square of the 

correlation coefficient approaches the maximum (step S47). 
25 Here, assuming that a film is formed on a semiconductor substrate (a fabrication 
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object surface) having a certain pattern, the height distribution measuring device 20 
measures a height distribution at the stage where this film is already polished by CMP. 
The measured result is measured data D22. On the other hand, the calculated 
two-dimensional distribution data of the irregularities on the fabrication object surface 
5 after CMP process is two-dimensional distribution data DP3. 

Therefore, for example, the polishing rate (A/sec) that is determined by the 
material of the polishing pad, the number of revolutions of the polishing pad and the like, 
the stress value (pascal) exerted on the polishing pad in the calculation object area, and 
the polishing time (sec), which were set in step S453, are changed so as to approach the 
10 measured data D22, i.e., increase the correlation coefficient. This is one of the 
parameter fittings after CMP process. 
C-2. Effects 

According to the method for simulating a CMP process and the simulator for a 
CMP process in the third preferred embodiment, before and after CMP process, a 
15 measured data and a simulation data are compared to obtain a correlation. Accordingly, 
the adjustment of parameters before CMP process can be separated clearly from the 
adjustment of parameters after CMP process, thus leading to the simulator that can 
consider various parameters. 

This provides the characteristic feature that if process conditions is changed or a 
20 new apparatus is added, adjustment may be accomplished merely by making a fine 
adjustment of parameters. 

In addition, various parameters such as the polishing rate, stress value exerted 
on the polishing pad and polishing time that are used for the calculation about polishing 
can be inspected by monitoring the correlation between the measured data after CMP 
25 process D22 and the calculated two-dimensional distribution data DP3 about the 
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irregularities on the fabrication object surface after CMP process. 
D. Fourth Preferred Embodiment 

Following is a method for simulating a CMP process and a simulator for a CMP 
process according to a fourth preferred embodiment of the present invention. Similar 
5 reference numerals have been used in the following figures to denote similar parts that 
have the same configuration as that shown in Figs. 1 and 2, and need not be described 
herein. 

D-l. Simulation Method and Simulator for CMP Process 

A method for simulating a CMP process will be described by referring to the 

10 configuration of a simulator for a CMP process 4 shown in Fig, 14, and by using the 
flowchart shown in Fig. 13. 

Referring to Fig. 14, the simulator 4 receives a pattern density data Dl per 
process from a pattern density data storage device 10, and receives a measured data D3 
about a height distribution of an under laid layer and height distributions before and after 

15 CMP per process from a height distribution measuring device 20. 

Based on a coordinate data contained in the pattern density data Dl, a pattern 
density two-dimensional distribution calculating part 41 1 expands the pattern density data 
such that a mesh data is arrayed in two dimensions to obtain a two-dimensional image. 
This provides a two-dimensional distribution image DP per process (step S51). 

20 In the measured data D3 about the height distribution of the irregularities on the 

under laid layer and the height distributions of the irregularities on the semiconductor 
device before and after performing a CMP process, which are provided from the height 
distribution measuring device 20, the condition of the irregularities on the semiconductor 
device is provided as a two-dimensional distribution image. However, the mesh of the 

25 pattern density two-dimensional distribution image DP is not always identical with that of 



31 

the height distribution measured data D3. In order to adjust such that the 
two-dimensional distribution image DP and measured data D3 have the same mesh, a 
mesh adjusting part 412 makes a mesh adjustment (step S52). 

Through the foregoing mesh adjustment operation, the mesh of the pattern 
5 density two-dimensional distribution image DP matches the mesh of the height 
distribution measured data D3, and these two data can be compared to each other. 

The mesh adjustment operation in the mesh adjusting part 412 is the same as 
that of the mesh adjusting part 112 shown in Fig. 2, and need not be described herein. 

At the completion of the mesh adjustment in step S 5 2, the resulting measured 
10 data D31 before it is subjected to a CMP process (before polishing), the measured data 
after CMP process (after polishing) D32 and the pattern density two-dimensional 
distribution image DP are provided to a height distribution calculating part 413. 

The height distribution calculating part 413 calculates a height distribution 
based on the pattern density two-dimensional distribution image DP, to obtain a height 
15 distribution data DPI about the fabrication object surface before CMP process (step S53). 
The method for calculating the height distribution based on the pattern density 
two-dimensional distribution image DP is described previously with reference to Figs. 5 
and 6, and need not be described herein. 

Subsequently, a measured data calculating part 414 adds data of a height 
20 distribution of an under laid layer to the height distribution data DPI about the fabrication 
object surface before CMP process, to obtain a height distribution data DPI 1 containing 
the height distribution data of the under laid layer (step S54). 

The height distribution data DPI 1 and measured data before CMP process D31 
are provided to a correlation coefficient calculating part 416, and these data and the 
25 measured data after CMP process D32 are provided to a CMP image calculating part 415. 



32 

Following is the operation of parameter fitting before CMP process. 
The correlation coefficient calculating part 416 performs a least squares 
analysis of the measured data D31 before CMP process and height distribution data DPI 1, 
to calculate a correlation coefficient (step S55). The analysis operation of step S55 is 
5 the same as that of step S4 shown in Fig. 1, and need not to be described herein. 

Here, assuming that a film is formed on a semiconductor substrate (a fabrication 
object surface) having a certain pattern, the thickness of the formed film is a fitting 
parameter. The height distribution measuring device 20 measures a height distribution 
at the stage where this film is formed on the semiconductor substrate. The measured 
10 result is measured data D31. On the other hand, the calculated value based on a pattern 
density data when forming this film is height distribution data DPI 1. 

Therefore, in a parameter fitting part 417 the formation thickness d2 of a 
laminated film SFM that was set in step S53 is changed so as to approach the measured 
data D31, i.e., increase the correlation coefficient (step S58). This is one of the 
15 parameter fittings before CMP process. 

Following is the operation of parameter fitting after CMP process. 
The CMP image calculating part 415 calculates a two-dimensional distribution 
data after polishing DP4 that is obtained by CMP process, from the height distribution 
data DP11 about the fabrication object surface before CMP process that contains the 
20 height distribution data of the under laid layer, as well as mechanical parameters such as 
Young's modulus and elastic coefficients of a polishing pad used in the CMP process 
(step S56). The operation of taking a CMP image in step S56 is the same as that 
described with reference to Figs. 1 1 and 12, and need not be described herein. 

The two-dimensional distribution data DP4 and measured data after CMP 
25 process D32 are then provided to the correlation coefficient calculating part 416. 
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The correlation coefficient calculating part 416 performs a least squares 
analysis of the measured data D32 and two-dimensional distribution data DP4, to 
calculate a correlation coefficient (step S57). The analysis operation of step S57 is the 
same as that of step S4 5 and need not to be described herein. 
5 Using the correlation coefficient obtained in step S57, as an index, the 

parameter fitting part 417 performs a parameter fitting such that the square of the 
correlation coefficient approaches the maximum (step S58). 

Here, assuming that a film is formed on a semiconductor substrate (a fabrication 
object surface) having a certain pattern, the height distribution measuring device 20 

10 measures a height distribution at the stage where this film is already polished by CMP. 
The measured result is measured data D32. On the other hand, the calculated 
two-dimensional distribution data of the irregularities on the fabrication object surface 
after CMP process is two-dimensional distribution data DP4. 

Therefore, for example, some parameters such as the polishing rate (A/sec) that 

15 is determined by the material of the polishing pad, the number of revolutions of the 
polishing pad, and the like, the stress value (pascal) exerted on the polishing pad in the 
calculation object area, and the polishing time (sec), which were set in step S453 (Fig. 1 1), 
are changed so as to approach the measured data D32, i.e., increase the correlation 
coefficient. This is one of the parameter fittings after CMP process. 

20 D-2. Effects 

According to the method for simulating a CMP process and the simulator for a 
CMP process in the fourth preferred embodiment, before and after CMP process, a 
measured data and a simulation data are compared to obtain a correlation. Accordingly, 
the adjustment of parameters before CMP process can be separated clearly from the 

25 adjustment of parameters after CMP process, thus leading to the simulator that can 
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consider various parameters. 

This provides the characteristic feature that if process conditions is changed or a 
new apparatus is added, adjustment may be accomplished merely by making a fine 
adjustment of parameters. 

In addition, various parameters such as the polishing rate, stress value exerted 
on the polishing pad and polishing time that are used for the calculation about polishing 
can be inspected by monitoring the correlation between the measured data after CMP 
process D22 and the calculated two-dimensional distribution data DP3 about the 
irregularities on the fabrication object surface after CMP process. 

Furthermore, with the configuration that the measured data of the height 
distribution of the under laid layer in the individual process is added to the pattern density 
data, it is possible to consider the influence of the previous process and therefore permit a 
simulation suitable for manufacturing semiconductor devices having a laminated 
structure. 

While the invention has been shown and described in detail, the foregoing 
description is in all aspects illustrative and not restrictive. It is therefore understood that 
numerous modifications and variations can be devised without departing from the scope 
of the invention. 



